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Abstract. We study the final state interaction in Stt decay of meson resonances at the Excited 
Baryon Analysis Center (EBAC) of JLab. We apply the dynamical coupled-channels formulation 
which has been extensively used by EBAC to extract N* information. The formulation satisfies the 
3k unitarity condition which has been missed in the existing works with the isobar models. We 
report the effect of the 37t unitarity on the meson resonance pole positions and Dalitz plot. 
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INTRODUCTION 

Why do we care the three-pion (or nnK, etc.) scattering system ? Some of recent and 
forthcoming experiments would call for serious study of the system to extract what they 
expect to find. One of such experiments aims to find the so-called exotic mesons in a 
reaction like ;r(or y)N M*N — )■ TtnnN [1, 2], where M* is an intermediate excited 
meson that could be exotic. An exotic meson lies outside of the constituent quark model 
and speculated to be a tetraquark state or hybrid state. The spectroscopy of the hybrid 
states provides information about gluon self-interactions. Another experiment seeks 
for physics beyond the Standard Model contributing to the CP violation in B- or D- 
decays [3, 4]. Both of the experiments involve heavy meson decays at very short-distance 
followed by the formation of lighter hadrons such as pions and/or kaons. And those light 
hadrons are to be detected after a great number of rescatterings. Thus, for extracting what 
has happened at the short-distance, it is essential to understand or precisely model the 
final state interactions of the light hadrons. 

What has conventionally been used for Dalitz plot analyses of experiments of this kind 
is the so-called isobar model in which it is assumed that two of the three mesons form 
an isobar state (/o,p, A'*, etc.), and the rest does not interact with the others (spectator). 
The three-body unitarity is obviously missing, and that is what we will address here. 
(Note that most isobar models do not take care of even two-body unitarity for the paired 
mesons.) Recently, the Excited Baryon Analysis Center of Jefferson Lab has extended 
their dynamical coupled-channels model [5], which has been extensively used to study 
baryon resonances, to three-meson scattering [6]. In this report, we describe our model 
and then apply it to the calculation of ;r2(1670) and ;r2(2100) pole positions and their 
371 decays (Dalitz plots). We will examine how significant unitarity is to extract the 
properties of these heavy mesons from the Dalitz plots. 
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FIGURE 1. K-R interaction: (a) 3k Z-graph; (b) M*-term 



COUPLED-CHANNELS MODEL 

We consider the decay of a heavy meson (M* -^3n) and a heavy meson resonance 
formed in 3n scattering (3;r — )■ M* 37r). Interesting quantities to calculate are the 
decay amplitude (and the corresponding Dalitz plot) for the decay, and the pole position 
of the resonance. In any case, we need 37t scattering amplitudes. For calculating 37t 
amplitude, first we need to develop a more basic kk model. In the following, we discuss 
our TtTt model and then calculation of the 3k amplitude. 

TtTt model For a simpler calculation, we employ an isobar-type model for TtTt interaction 
so that TtTt ^ R ^ TtTt (or KK), where i? (= /o,p or /2) is an isobar. The TtTt potentials 
for a partial wave (orbital angular momentum L, total isospin /) are parametrized as 
T.RfRi{p'){W —mR)^^f]^^j{p), where W and mr are the total energy and the bare mass 
of the isobar R, respectively; fR,i{p) is the R -H^ / (/ = TtTt or KK) vertex with p being 
the relative momentum of TtTt (KK). The TtTt scattering amplitude is obtained by solving 
the coupled-channels Lippmann-Schwinger equation with this potential. 

;r-isobar scattering equation Because we have introduced the isobars (R), the partial 
wave amplitude specified by (total spin and parity) and T (total isospin) can be 
obtained by solving the corresponding Tt—R scattering equation: 



Tip{p',P;W)=ViJip',p;W)+Y, / q^dqViJ{p',q;W) GyY{q,W)T^J {q,p;W) , (1) 
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where the indices «,/3,7 specify channels (tt/q or np or 71/2), Ey{q) = A/m^. + ^2_j_ 
^m^^ + q'^. The Tt — R Green function, Gyyiq.W), is given by [G^^(^, V7)]yy = [V7 — 
£'y(^)]5yy' — Eyy W), with the self-energy (Eyy) defined as 

i V ^^r(P)E^Y iP) ^Ef{q) +p^W- E„{p) - ^Ef{q) +p^ + ie 

The potential, V^''p(p',p;W), includes Z-graphs [Fig. 1(a)] derived from the TtTt model 

developed above. (We ignore Z-graphs with TtKK intermediate states for simplicity.) 

The Z-graphs, together with the self-energy [Eq. (2)], are essential to maintain the 3Tt 
unitarity. One may also include M*-term [Fig. 1(b)] for the K-R potential. 




FIGURE 2. M*-decay amplitude: The bulb labeled 't' is the T-matrix calculated with Eq. (1) in which 
the potential includes only the Z-graphs. The gray circle indicates the self-energy of the isobar [Eq. (2)]. 



APPLICATIONS AND RESULTS 

Effect of Z-graph on pole position We examine the effect of Z-graph on the pole 
position of ;r2(1670) and ;r2(2100). The procedure for finding the pole position suitable 
for our model has been developed [7]. A pole position is at a complex energy (W) for 
which W — M^f* - ^M* (W) = is satisfied; M%* is the M* bare mass and Lm* (w) is the 
self energy of M* given by 

POO 

Em*(W) = J^/ q^dq gM*,a{^)Gaii{q,W)gM*,p{q) , (3) 

where gM*,ji (q) is a M* — )■ kR vertex, and gM*,a{q) is a dressed vertex defined as 

f°° p 
8M*,a{p) = 8M*,a{p)+Y. / ^^^^ 8M*,y{q) Gyy{W,q) t^yJ^{q,p;W) , (4) 

fpf 

where ti „ is the T-matrix calculated with Eq. (1) in which the potential includes only 

the Z-graphs. We fit M^* and gM*,a to pole positions and branching ratios of K2's; data 
are taken from the Particle Data Group. Then we eliminate the Z-graphs from our model, 
that is, we calculate Em*(VF) with Eq. (3) in which the dressed vertex is replaced by the 
bare one. We find the pole positions again from W — M^* — Zm*{W) = 0. The pole for 
;r2(1670) changes from 1672 - 130/ MeV to 1689 - 133/ MeV, and for ;r2(2100), from 
2090 - 313/ MeV to 2084 - 346/ MeV. Thus, the Z-graphs change the pole position by 
~ 10% [2 ~ 3%] for ;r2(2100) [;r2(1670)]. 

Effect of Z-graph on Dalitz plot We calculate the decay amplitude of a M* as graphi- 
cally shown in Fig. 2, and then obtain the corresponding Dalitz plot. In the first term of 
Fig. 2, two pions are paired, and the 3rd pion is the spectator. The second term includes 
rescatterings due to the Z-graphs. Here, we show the Dalitz plot for the ;r2(2100) decay 
in Fig. 3 (left). In order to examine the effect of the Z-graphs, we turn off the Z-graphs, 
i.e., we calculate only the first term in Fig. 2. We show the ratio of the Dalitz plots with 
and without the Z-graphs in Fig. 3 (right). In general, the Z-graphs change both the mag- 
nitude and the shape of the Dalitz plot. The effect is more significantly seen in the Dalitz 
plot than in the shift of the pole position. 

We may fit to the Dalitz plot of Fig. 3 (left) with our model without the Z-graphs, by 
varying the M* — )> ttR couplings and cutoff. With a tentative error of 5% assigned to the 
Dalitz plot, we achieved the fit with x^/ data ~ 1 . 1 . It would be interesting to examine the 



TABLE 1. The properties of ;r2(2100) extracted from models with (2nd row) and without (3rd row) 
the Z-graphs. The second column is the width at the pole, 3-5th columns are the branching ratios. 
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FIGURE 3. (left) Dalitz plot of (2100) decay; (right) Ratio of DaUtz plots with and w/o the Z-graphs. 



properties of M* extracted from the two models (with or without the Z-graphs) both of 
which reproduce the same Dalitz plot. In Table 1, we show the width and the branching 
ratios of ;r2(2100) extracted from the two models. We find a significant difference in 
the properties from the two models. Although the effect of the Z-graphs depends on a 
system, still this results indicates an importance of including the Z-graph (371 unitarity) 
in analyzing the Dalitz plot. 
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